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ABSTRACT 

We perform multi-dimensional radiative transfer simulations to compute spectra for a hydro- 
dynamical simulation of a line-driven accretion disk wind from an active galactic nucleus. 
The synthetic spectra confirm expectations from parameterized models that a disk wind can 
imprint a wide variety of spectroscopic signatures including narrow absorption lines, broad 
emission lines and a Compton hump. The formation of these features is complex with con- 
tributions originating from many of the different structures present in the hydrodynamical 
simulation. In particular, spectral features are shaped both by gas in a successfully launched 
outflow and in complex flows where material is lifted out of the disk plane but ultimately 
falls back. We also confirm that the strong Fe Ka line can develop a weak, red-skewed line 
wing as a result of Compton scattering in the outflow. In addition, we demonstrate that X-ray 
radiation scattered and reprocessed in the flow has a pivotal part in both the spectrum forma- 
tion and determining the ionization conditions in the wind. We find that scattered radiation 
is rather effective in ionizing gas which is shielded from direct irradiation from the central 
source. This effect likely makes the successful launching of a massive disk wind somewhat 
more challenging and should be considered in future wind simulations. 

Key words: radiative transfer - methods: numerical - galaxies: active - X-rays: galaxies 



1 INTRODUCTION to the escape speed from the stellar photosphere (e .g. ICastor et al.l 

ll975l ; lFriend & Abbott|[T98rj: IPauldrach et alJI 1986b . For AGN ac- 
Estimates for the black hole mass (M bh ) and bolometric luminosity cretion disks, launching of a wind is hindered by the large lumi- 
(Lboi) of the most luminous active galactic nuclei (AGN) suggest nos j ty f ionizing X-ray radiation produced in the vicinity of the 
that they are accreting close to the Eddington rate, M Edd . Under central black hole. Line-driving becomes rapi dly less efficient as 
such conditions, radiation forces should play a role in determin- the ionization state of t he gas is raised (see e.g. lStevens & Kallmanl 
ing the properties of the accretion disk atmosphere. In particular, |l99d : lArav et aljl994l) . As a consequence, for luminosities around 
since much of the luminosity of the inner regions of the disk is rj.OOl < L/L Ed d ^ 0.1, line-driving is unable to launch a success- 
radiated in the ultraviolet (uv) part of the spectrum, a significant f u j wmcl but should give rise to a topological^ complex "failed 
radiative force should arise owing to the large number of line tran- wind " m wnich mater ial is lifted out of the disk plane but, due to 
sitions which occur in this spectral region. By analogy to hot star over-ionization by t he X- ray radiation, does not reach e scape speed 
atmospheres, it is to be expected that radiation pressure due to ul- jProga & Kallmanl |2004 hereafter PK04, see also IProga et al.l 
traviolet line driving can e xceed that due to e l ectron scattering by |200Ct|Progj2005b . Above L ~ 0.1 L Edd , however, numerical sim- 
factors up to M m ^ ~ 1 3 dCastor et al.||l975|; |Stevens & Kallmanl ulations suggest that even the ionizing X-ray flux may be unable to 
1 1990]; |Proga et al.||2003) under optimal conditions. Thus the radi- f u n y suppress the launching of a line-driven disk wind. In particu- 
ation force due to spectral lines might be expected to affect the PK 04 performed radiation-hydrodynamical simulations for an 
structure of a uv-bright accretion disk for luminosities as small as accretion disk around a supermassive black hole (M bh = 10 s M ). 
~ Af max LEdd- They investigated the case of L ~ 0.5L E dd assuming that 10 per 
In luminous hot stars line-driving by the uv radiation field cent of the luminosity was radiated as X-rays from a centrally con- 
is able to launch fast winds with terminal velocities comparable centrated source. This central source could ionize the gas but was 
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not included as a source of radiation pressure - the line driving 
force was supplied by the uv radiation of the accretion disk. Their 
simulations predicted that, despite the ionizing effects of the X-ray 
radiation, a significant mass flux was successfully accelerated that 
formed an outflow with terminal velocity comparable to the escape 
speed from the inner accretion disk. 

The flow structure found by PK04 is complex and significantly 
time- variable but it is persistent and characterized by three regimes: 
a cold dense equatorial outflow, a hot diffuse polar flow and a thin 
layer of moderate density rapidly outflowing warm material at in- 
termediate inclination. The typical mass-loss rate of the flow was 
~ 0.1 — 0.2 Mq yr" 1 , this being dominated by the warm, rapid out- 
flow component. In addition, the simulation predicts a dense "failed 
wind" region above the inner disk which plays a pivotal part in de- 
termining the wind structure owing to both the ram pressure it ex- 
erts and its potential role in shielding the gas in the outer regions 
of the flow from the ionizing X-ray radiation. We note, however, 
that PK04 assumed that all the X-ray radiation is emitted by a point 
source and they neglected scattering or reprocessing of X-rays in 
the flow. As will be discussed later, this likely overestimates the 
ease with which the "failed wind" can shield the outer parts of the 
flow. 

Given the geometrical and kinematic complexity of the PK04 
wind structure, predicting its spectroscopic signatures is challeng- 
ing. But given the large mass-loss rate of the warm stream and the 
significant covering fraction of the "failed wind", it is certainly ex- 
pected that observational signatures will arise. Since a large frac- 
tion of the material around the inner parts of the flow should be 
strongly ionized, the clearest signatures of the absorption and emis- 
sion from t he wind structure ma y be expected to manifest in the 
X-ray band. lSchurch et al] 12009) performed radiation transfer cal- 
culations for snapshots of the wind structure obtained by PK04, 
again assuming that all X-rays are emitted by a single point source. 
They confirmed that, depending on the observer's orientation, the 
wind structures could strongly absorb X-ray radiation and there- 
fore imprint observab le signatures. Although an important step, the 
Schur ch et al. I d2009l) study suffers fro m the shortcoming that the 
radiative transfer code used (XSCORT; ISchurch & Pond 12007b is 
only one-dimensional. Thus, although they could study the absorp- 
tion properties of the disk wind, they did not realistically address 
the role of the wind in producing/shaping emission features. 

We have developed a multi-dimensional Monte Carlo radia- 
tive transfer code which is able to compute orientation-dependent 
spectra for wind models taking proper account of s cattering and re- 
processing of ra diation in the ou tflow ( Sim 2005', ISim et ai]|2008l ; 
hereafter Paper I. ISim et alj20 10; hereafter Paper II). Applying this 
code to simply-parametrized disk wind geometries has shown that 
scattered/reprocessed radiation has a very important role in the for- 
mation of the X-ray spectrum and that outflow models may have 
the capacity to explain a wide variety of X-ray spectroscopic fea- 
tures (see Papers I and II). In particular, highly blue-shifted ab- 
sorption features, as have been reported in a num ber of luminous 
AGN (see section 6.1 of iTurner & Miller! 120091 for a recent re- 
view), appear for a subset of inclination angles while Fe Ka emis- 
sion is predicted for all observer orientations. The line emission is 
typically broad and can develop somewhat red-skewed line wings 
(Papers I and II) owing to the effects of Compton scattering in 
the flow (see e.g. iLaming & Tita rchuk 20041; [Laurent & Titarchukl 
l2007l ; lTitarchuk et al.ll2009l) . 

To date, however, we have only applied our Monte Carlo code 
to parametrized outflow models with imposed mass-loss rates and 
flow geometries. Although such simple models provide a conve- 



nient means to investigate the formation of the X-ray spectrum in 
an outflow, they do not adequately represent the full complexity 
of the wind structures found in the radiation-hydrodynamical sim- 
ulations mentioned above. In particular, they adopt a very simple 
velocity law (which is not dynamically consistent with any par- 
ticular launching/acceleration mechanism) and assume that all the 
material that is lifted out of the disk plane is part of a flow that suc- 
cessfully escapes from the system. Thus, at best, the parametrized 
wind models might represent components similar to the warm out- 
flow in the PK04 simulation but they neglect any other structures 
which accompany it (e.g. the "failed wind"). It is therefore impor- 
tant to investigate whether there are characteristic differences be- 
tween the X-ray spectra obtained from simplified wind models and 
dynamically self-consistent flow models since these are relevant to 
the quantitative interpretation of observations. 

Here we take the step of connecting our studies of the forma- 
tion of X-ray spectra in multi-dimensional flows (Papers I and II) to 
a physically motivated disk wind model (PK04). The objectives of 
this study are two fold. First, by computing synthetic X-ray spectra 
for a hydrodynamical simulation we can more reliably investigate 
which spectral features can plausibly form in disk winds and as- 
sess how well synthetic spectra from parametrized models (such as 
those used in Paper II) correspond to what might be expected from 
more realistic flow structures. Second, by using a radiative trans- 
fer scheme that readily incorporated true absorption, scattering and 
reprocessing of radiation in the wind, we can begin to assess the im- 
portance of scattered radiation for the physical processes at work in 
a line-driven wind. We begin, in Section[2]with a description of the 
setup of our simulations. The results are presented in Section[3]and 
their implications discussed in Section|4] 



2 MODEL SETUP 

2.1 Radiation-hydrodynamical wind simulation 

To provide a set of self-consistent gas properties for our Monte 
Carlo radiative transfer simulations, we use the hydrodynamic ac- 
cretion disk wind simulations presented by PK04. The hydrody- 
namical simulations assume a 10 8 Mq black hole accreting at 
1.8 Mq yr" 1 . These parameters correspond to an accretion rate 
of 0.5 AfEdd assuming a radiative efficiency of 6 per cent. The gas 
is assumed to accrete through an optically thick, geometrically thin 
standard disk. The disk thermal uv radiation drives the wind. The 
simulations also include Comptonized continuum radiation from a 
central source which represents the emission from the inner parts 
of the accretion disk that are not covered by the computational do- 
main. This component is not included in the line-driving calcula- 
tion for the wind, because in the X-ray band there are relatively 
few line transitions. However, the central X-ray component has a 
major effect on the derived photoionization state. The properties of 
the outflow are calculated over the radial range 60 - 3000 r g (where 
r g — GMbh/c 2 is the gravitational radius of the black hole). 

To solve the equations of gas dynamics, PK04 used the ZEUS 
code which has been tested extensively (Stone & Norman 1992). 
Their version of the code has been modified to implement the radi- 
ation force, and radiative heating and cooling. PK04 approximated 
the radiative accelerat ion due to lines u sing a modification of the 
method developed bv lCastoretd]jl975h . The radiative driving due 
to line scattering is treated using the generalized Sobolev approxi- 
mation. This assumes the transport in the line core is dominated by 
the local velocity gradient. The coupling between the dynamics and 



3 



the radiation field is provided by radiation pressure and by thermal 
pressure due to radiative heating. The radiation pressure depends 
on the gas opacity, which in turn depends on the ionization balance. 
The effects of X-ray ionization on the radiati on force due to lines 
is calc ulated using the procedure outlined by IStevens & Kallman 
1 1990). In this procedure, the "force multiplier" of ICastor et al 
(1975) (the ratio of the full radiation force to that due solely to 
electron scattering) is suppressed by a factor which depends on the 
ionization parameter and temperature. The effects of heating by ra- 
diation are also taken into account by these calculations. Specifi- 
cally, PK04 calculated the gas temperature assuming that the gas 
is optically thin to its own cooling radiation. Thus the net cooling 
rate depends on the density, the temperature, the ionization param- 
eter, and the characteristic temperature of the X-ray radiation. In 
this case it is possible to fit analytical formulae to the heating and 
cooling rate obtained from detailed photoionization calculations for 
various gas parame ters. PK04 used a fit to photoionization calcu- 
lations obtained by Blondin ( 1994) who included Compton heat- 
ing/cooling, X-ray photoionization heating/recombination cooling, 
bremsstrahlung and line cooling. 

2.2 Monte Carlo radiative transfer simulations 

The hydrodynamical simulations described above provide the time- 
dependent density and velocity structure of the disk wind. To com- 
pute synthetic spectra from these, we have performed radiative 
transfer simulations using the Monte Carlo code described in Pa- 
pers I and II. 

For the radiative transfer simulations, the Monte Carlo code 
was modified to accept a generalized axisymmetric wind with den- 
sity and velocity specified via input data. As for the parametrized 
wind models adopted in Papers I and II, these input data define a 
2D grid of wind properties (in this case, with properties depend- 
ing on polar coordinates r and 8). The mass density in each wind 
grid cell is assumed to be uniform and taken from the hydrody- 
namical model. The three components of velocity (v r , vg and v^) 
at every point in the wind are obtained by linear interpolation (in r 
and 9) between the values which are provided at the boundaries of 
each wind grid cell. This interpolation is necessary since use of the 
Sobolev approximation for line transitions requires that the velocity 
is everywhere a smooth function. 

In certain regions, particularly close to the disk plane where 
the transition between the disk atmosphere and the outflow lies, 
the hydrodynamical simulations predict rather high densities. In 
some cases, this makes certain grid cells very optically thick. Phys- 
ically, any X-ray photons which penetrate deep within these opti- 
cally thick layers will be thermalized and contribute to heating of 
the disk atmosphere. In the parametrized models considered in Pa- 
per II, this sink of X-ray photons was roughly accounted for by 
assuming that all Monte Carlo quanta which reach the xy-plane 
in the simulation strike the optically thick disk and are then lost. 
That approach neglects reflection by the disk but avoids the need to 
track the Monte Carlo quanta as they propagate through very op- 
tically thick material. In the more realistic model considered here, 
however, there is no sharply-defined boundary between the disk at- 
mosphere and the wind. Thus reflection by the dense material at 
the very base of the outflow is automatically included in the sim- 
ulations. However, this introduces additional computational chal- 
lenges. In particular, quanta occasionally propagate deep into the 
optically thick regions where they interact many times and have 
very low probability of reemerging without being thermalized (and 
therefore lost to the X-ray regime since the expected temperature of 



the disk atmosphere is relatively low). To deal with these quanta, a 
cut is introduced whereby the flight paths of Monte Carlo quanta 
which propagate sufficiently deep into the base of the wind are 
terminated and it is assumed that these quanta contribute nothing 
further to the X-ray spectra or ionization state of the wind. A very 
conservative cut has been applied in the simulations presented here: 
quanta are terminated if they penetrate deep enough into the base 
of the wind that the Compton optical depth to re-emerge from the 
wind is greater than twenty in all directions. 

As described in Paper II, the ionization and thermal structure 
of the wind are determined iteratively from the radiation field prop- 
erties and the assumptions of ionization and thermal equilibrium. 
The ionization balance accounts for photoionization (including K- 
/L-shell ionization followed by ejection of Auger electrons), colli- 
sional ionization, radiative recombination and di-electronic recom- 
bination. Heating by photoionization, Compton down-scattering 
of X-ray photons and free-free absorption are balanced against 
bremsstrahlung, Compton cooling by low energy disk photons, 
bound-free recombination, bound-bound line transitions and cool- 
ing due to adiabatic expansion to estimate the local tempera- 
ture. Note that our treatment of Compton cooling remains par- 
ticularly approximate since we do not simulate the transport of 
the accretion disk photons in detail (see Paper II). The assump- 
tions of ionization/thermal equilibrium are not perfectly valid in 
regions of the wind where the density is sufficiently low that the 
recombination/cooling timescales are long compared to the flow 
timescale. However, it alleviates the need to introduce explicit time- 
dependence in the calculation - a major computational saving. In 
general, we would expect that the most important consequence of 
departures from equilibrium is that the ionization state becomes 
frozen-in in the outer, low-density flow regions since recombina- 
tion may become too slow to maintain local ionization equilibrium. 
As we shall show below, however, the equilibrium assumption al- 
ready leads to near fully-ionized conditions in much of the outer 
wind - thus accounting for any further effective reduction to the 
recombination rate is unlikely to dramatically alter the typical ion- 
ization state. 

The radiative transfer simulations were performed using the 
same set of atomic data as described in Paper II. This includes the 
K-shell ions of the astrophysically abundant metals, the L-shell ions 
of the important second and third row elements and the highest few 
M-shell ions of Fe and Ni. Note that this data set does not include 
M-shell ions of iron below Fe X - thus when the computed ioniza- 
tion state favours this ion it is likely that the true degree of ioniza- 
tion is lower. Similarly, a lower limit of log T e [K] =4.0 is imposed 
in the radiative transfer simulations since (i) the atomic data used 
is likely inadequate to describe the cooling to lower temperatures 
and (ii) the uv radiation of the disk - which is neglected here - 
should prevent the inner regions of the outflow from dropping to 
significantly lower temperatures. 

2.3 Simulation parameters 

The relevant physical and numerical parameters adopted in the 
simulations are given in Table [TJ Most of these are carried over 
or derived from the parameters adopte d in the radiation hydrody- 
namics simulation dProgaetal.ll2000L PK04). For this study we 
have selected two snapshots of the wind structure; specifically we 
chose the wind conditions from timesteps 800 and 955 from the 
PK04 simulation. These timesteps were chosen since they are well- 
separated in simulation time (timestep 955 is later than timestep 
800 by At ~ 5 years). This allows us to study two independent 
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Table 1. Parameters for the simulation. 



Parameter Value 

mass of central object, 10 8 Mq 

source luminosity, L(, i 6.5 X 10 45 ergs s —1 (~ 0.5Lridd) 

X-ray source luminosity (2 - 10 keV), L x 2.7 X 10 44 ergs s" 1 (~ 0.02L Edd ) 

source power-law photon index, r 2.1 

range of source photon energies in simulation 0.1-511 keV 

size of primary emission region, r er 6r 9 = 8.8 X 10 13 cm 

inner radius of disk.r^ 6r 9 = 8.8 X 10 13 cm 

inner radius of simulated wind region, r m j„ 60 r g = 8.8 X 10 13 cm 
outer radius of simulated wind region, Tmax 

3000 r„ 

outer radius of Monte Carlo Rad. Trans, simulation grid 5 X 10 16 cm 

3D Cartesian Rad. Trans, grid cells 256 X 256 X 256 

2D radiation-hydrodynamics polar grid zones 100 X 140 



realizations of the flow pattern and investigate how the spectral fea- 
tures should vary on timescales comparable to the typical interval 
between repeat observations of well-known AGN. Both snapshots 
show the characteristic structure of the PK04 simulation mentioned 
in Section[T](a hot polar flow, warm rapid outflow, cold equatorial 
flow and a "failed wind" region). The relative importance of these 
structures is slightly different between the two-snapshots, however: 
this can be seen in the upper right panels of Figure [T] (mass density 
of the two-snapshots) and in Figure[2] which compares the hydro- 
gen column densities for the two snapshots for different inclination 
angles, 6 (measured relative to the rotation (z-)axis.) 

For our radiative transfer simulations we need to fully specify 
the X-ray source radiation. We assumed that the radial extent of the 
primary emitting region is comparable to the size of the innermost 
stable orbit around a Schwarzschild black hole (r em = 6r 9 ) and 
we also assumed that a geometrically thin/optically thick accretion 
disk lies in the mid-plane of the simulation, extending from an inner 
radius of rd = 6r 9 to the outer boundary of the simulation grid. 
The X-ray source is assumed to radiate isotropically. 

We assume that the primary X-ray radiation source is a pure 
power-law across the photon energy interval we simulate (0.1 < 
hv < 511 keV; see Paper II). We have adopted a photon index 
of r = 2.1. This is reasonably typical for the class of objects for 
which disk winds have been discussed as a possible site for the for- 
mation of X-ray s pectral features (e.g. T ~ 1.8 for PG1211+143 
JPounds etal.ll2003ll . T ~ 2.4 for Markarian 766 iMiller et all 
l2007h or T ~ 2.2 for MCG-6-30-15 I Mil ler et alj|2008rl ). To es- 



timate the appropriate normalisation of the primary power-law X- 
ray emission, we used the source X-ray luminosity of 0.05 Lmi 
adopted by PK04. The role of the X-ray luminosity in the radiation 
hydrodynamics simulation is to estimate the ionization parameter 
required to compute the line driving force multiplier. PK04 used 
the relatio nship between force multip lier and ionization parameter 
derived bv lStevens & Kallmad (l990) which was based on calcula- 
tions assuming a bremsstrahlung source spectrum with temperature 
Tx ~ 10 keV. Therefore, we fixed the normalisation of our input 
X-ray luminosity spectrum (L(E) oc E~ F+1 ) to give the same to- 
tal X-ray flux (integrated from 1 to 20 keV) as for a bremsstrahlung 
spectrum (L(E) = 0.05L E dd exp(-E/T x )/kT x ). This yields a 2 
- 10 keV luminosity of Lx = 0.02Z/Edd which is comparable to 
the observationally motivated estimates for Lx /Leu adopted for 
the models described in Paper II. 

The Monte Carlo code requires that the total photon lumi- 
nosity (Lboi) be specified in order to estimate the role of Comp- 



ton cooling in the wind (see Paper II). We have adopted Lboi = 
6.5 X 10 45 ergs s _1 (= 0.5I/Edd) which is consistent with the ac- 
cretion rate (0.5AfEdd) and radiative efficiency (77 = 0.06) adopted 
by PK04. 



3 RESULTS 

As described in Paper II, the ionization state and kinetic tempera- 
ture in the wind are computed from the radiation field properties 
via an iterative sequence of Monte Carlo simulations. We will first 
describe the results of this part of the calculation and then present 
the synthetic spectra. 



3.1 Temperature structure 

Figure [T] shows the temperature distribution computed from the 
Monte Carlo radiative transfer simulations for the two snapshots 
(top and second row, left panels). In both cases the temperature 
structure is complex but can be easily understood as a consequence 
of the density distributions (top and second row, right panels in 
Figure [TJ and the central location of the primary X-ray source. 
The low-density polar wind cone is strongly irradiated by the X- 
ray source and is therefore both hot and very highly ionized (see 
Figure [3] discussed below). Below this region lies a warm stream 
(T e ~ 10 6 K) where the relatively dense portions of the wind are 
exposed to the primary X-rays. Deeper inside, the wind becomes 
significantly cooler, down to T e < 2 xlO 5 K in the outer regions. 
The densest material (comprising the "failed wind" region and the 
atmosphere of the outer parts of the disk) is not significantly heated 
by the X-ray source and remain at the lower temperature boundary 
imposed in the simulations (T e — 10 4 K). 

The computed temperature (and ionization) structure of the 
disk wind are noticeably different from those obtained in the origi- 
nal hydrodynamical simulations of PK04. For comparison, the tem- 
perature distribution for the timestep 955 snapshot obtained from 
those simulations is also shown in Figure Q] (bottom left panel). 
There are several qualitative differences: most importantly, the 
dense regions of the outflowing gas are significantly cooler that 
obtained here. This can be attributed to the differences in the radia- 
tive transfer schemes employed. In particular, PK04 did not include 
scattering - they assumed that the wind is a pure absorbing struc- 
ture. Our Monte Carlo method is considerably more sophisticated 
and includes both scattering and true absorption. To illustrate the 
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Figure 1. Distribution of kinetic temperature (T e ) computed from our Monte Carlo radiative transfer simulations (top/middle left) and mass density 
(top/middle right) for the two snapshots we consider from the PK04 simulation (timesteps 800 and 955; top and middle rows, respectively. For compari- 
son, we also show the temperature obtained for timestep 955 in the PK04 radiation-hydrodynamics simulations (lower left). The lower right panel shows the 
temperature distribution obtained with our Monte Carlo code (for the timestep 955 snapshot) when heating by scattered/reprocessed radiation is neglected. Note 
the logarithmic axes. The figures show the complete computational domain of the hydrodynamical simulation (60r 9 < r < 3000r 9 ). A lower temperature 
limit of T e = 10 4 K was imposed in the Monte Carlo radiative transfer simulations. 
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Figure 2. Hydrogen column density, JVjj versus inclination angle 6 (mea- 
sured relative to the rotation (z)-axis) for the two snapshots from the PK04 
simulation for which we compute spectra (filled circles for timestep 800, 
open circles for timestep 955). 

importance of this, we also computed the temperature distribution 
if only the heating due to direct irradiation is considered for the 
same snapshot (i.e. heating by radiation which is either scattered 
or reprocessed in the wind is neglected; see bottom right panel in 
Figure \Q. This reproduces the low temperatures found in the en- 
tire region shielded by the "failed wind" in the radiation hydro- 
dynamics simulation. Thus our Monte Carlo simulations show that 
scattered/reprocessed X-ray radiation can be rather important in de- 
termining the physical conditions in the disk wind, particularly for 
material that is shielded from direct radiation by the X-ray source. 

There are several other, relatively inconsequential, differences 
between the temperature distribution obtained from the Monte 
Carlo radiative transfer simulations and the radiation hydrodynam- 
ics. The polar wind is generally found to be hotter with the radiation 
hydrodynamics code than the Monte Carlo code. This is partly a 
consequence of the different forms of the adopted X-ray spectrum. 
In both codes, much of this material is predicted to be close to the 
Compton temperature T c but PK04 adopted an X-ray temperature 
of 10 keV (yielding T c ~ 3 x 10 7 K) while the input spectrum of 
the Monte Carlo simulation (power-law X-ray spectrum combined 
with total luminosity Lboi) gives a lower value of T c ~ 10 7 K. Also, 
in the radiation hydrodynamics calculation, the region above the 
dense disk wind contains structured regions in which the tempera- 
ture is much higher due to shock heating (interactions between ex- 
panding and falling parts of the flow). Since non-radiative heating 
is not included in our Monte Carlo code, such high temperatures are 
not found in our radiative transfer simulations. We note, however, 
that the exactly physical conditions in these polar regions are rather 
unimportant for the spectrum since this material is generally opti- 
cally thin. Thus the consequences of any discrepancy between the 
radiation-hydrodynamics and Monte Carlo radiative transfer simu- 
lations in the low-density polar regions are relatively minor. 

3.2 Ionization conditions and Fe Ka forming regions 

Since the X-rays are responsible for both ionization and heating 
in our Monte Carlo simulations, the distribution of ionization state 
in the wind is qualitatively similar to the temperature distribution. 



This is illustrated in Figure|3]which shows the distribution of mean 
iron ionization state computed for both snapshots. In the low den- 
sity polar regions the material is almost fully ionized whilst the 
denser streams in the flow are less ionized. Only the dense outer 
regions of the disk atmosphere (large radial coordinate) and the 
"failed wind" region are not significantly ionized. 

A noteworthy convenience of the Monte Carlo method is that, 
by analyzing the flight paths of quanta, it is relatively easy to iden- 
tify which parts of the outflow are responsible for the formation 
of particular spectral features. Since much of the discussion below 
will be devoted to the Fe K region, we have examined the Monte 
Carlo quanta which contribute to this part of the spectrum in detail. 
Specifically, since the Monte Carlo simulation follows the com- 
plete flight paths of the quanta, we can flag each quantum to iden- 
tify which physical processes (e.g. bound-bound transition, Comp- 
ton scattering etc.) were responsible for determining its properties 
(e.g. photon energy, direction of emission etc.) as it left the com- 
putational domain. To understand the formation of the Fe K band 
spectra, we identified all escaping quanta whose last physical inter- 
action (excluding Compton scattering) resulted in line emission via 
an Fe Ka transition. These packets are responsible for Fe Ka emis- 
sion features in the spectrum. We then binned these quanta based 
on the position of the Fe Ka emission which gave rise to them. The 
resulting distribution, which is shown for the PK04 timestep 800 
snapshot in Figure [4] indicates which portions of the flow are re- 
sponsible for the creation of the Fe Ka emission features. Since 
the emission line profile can be affected by Compton scattering, we 
also recorded the positions of the final Compton scattering event for 
the quanta flagged as contributing to the Fe Ka emission. These al- 
low us to identify the parts of the model in which scattering of the 
Fe Ka emission takes place - this region is also indicated in Fig- 
ure [4] Note that the plot does not retain any information related 
to observer orientation - all escaping Monte Carlo quanta are in- 
cluded, independent of their direction of propagation. 

Much of the Fe Ka line formation is concentrated in the in- 
ner parts of the flow structure: roughly 50 per cent of the escaping 
Fe Ka photons are created within a radial extent of r < 2 x 10 15 cm 
while ~ 90 per cent are formed inside r < 7 x 10 15 cm. Compar- 
ing Figure [4] to Figures Q] and [3] it is apparent that the most in- 
tense region of Fe Ka formation is around the inner surface of the 
"failed wind". This is expected since it is here that the X-ray radi- 
ation from the central source most directly strikes relatively dense 
material from which it can be reflected. Ka formation continues 
throughout most of the moderate density outflowing region outside 
the "failed wind", however, and there is also a reflection component 
formed from the dense material near the disk atmosphere (equato- 
rial regions) at large radii. Shielding by the "failed wind" means 
that this region is not directly irradiated by the central source but 
reflection still occurs thanks to illumination by X-rays scattered in 
the upper portions of the wind. 

In the outer wind, there are regions where the ionization state 
is too high for significant numbers of Fe Ka photons to be pro- 
duced but where Compton scattering of Ka photons created deeper 
in the wind can still occur. Regions such as this are potentially crit- 
ical since Compton scattering in a fas t outflow may have a role i n 
shaping emission line profiles (see e.g. lLaurent & Titarchukll2007h . 
In our simulations, escaping Fe Ka photon packets undergo an av- 
erage of about two Compton scatterings between their creation and 
escape (~ 40 per cent escape with no Compton scattering while 
~ 15 per cent undergo five or more Compton scattering events). 



7 





10 15 10 16 
x (cm) 



Figure 4. Regions responsible for the creation and subsequent scattering 
of the Fe Ka photons which contribute to the spectra for the snapshot from 
timestep 800 of the PK04 simulation. Each grid cell is shaded orange/red if 
Monte Carlo quanta representing Fe Ka photons were created in that cell 
and subsequently escaped without undergoing any further physical inter- 
actions (except Compton scattering). The intensity of the shading is pro- 
portional to the Fe Ka luminosity of the cell (ergs s — 1 ; i.e. the volumes 
of the cells are not taken into account). Grid cells are shaded blue if no 
Monte Carlo quanta representing Fe Ka photons were created in that cell 
but Compton scattering of escaping Fe Ka photons did occur there. White 
regions indicate that either no Fe Ka photons were created/scattered there 
or that any such photons were subsequently destroyed and so did not con- 
tribute to the spectra. 



3.3 X-ray spectra 

3.3.1 General characteristics 

Given the complex structure of the outflow (see FigureQ} and the 
high column densities for many lines-of-sight (see Figure |2j the 



emergent spectrum is expected to be strongly dependent on the ob- 
server inclination angle (8 = cos" 1 /i). As in Paper II, we extracted 
spectra for twenty different observer orientations, uniformly sam- 
pling the interval 1 > fi > 0. Figure [5] shows the spectra com- 
puted for three example observer orientations for the snapshot from 
timestep 800 of the PK04 simulation. Despite the underlying com- 
plexity of the wind model, these spectra are qualitatively rather sim- 
ilar to those obtained from the simply-parameterized smooth wind 
models used in Paper II. The spectra can be broadly classified into 
three groups, differentiated by the observer orientation; these are 
exemplified by the spectra shown in Figure|5]and discussed below. 

For relatively pole-on inclinations (first panel of Figure|5]l, the 
direct component of radiation is largely unattenuated and supple- 
mented by reflected radiation from the wind. As described in Sec- 
tion 13.21 (see Figure [4), the most intense reflection occurs at the 
inner irradiated surface of the wind (the "failed wind" region) and 
this gives rise to a reflection component when observed pole on. 
Qualitatively similar to reflection contributions in the parametrized 
disk wind models us ed in Paper II and disk reflection models 
dRoss & Fabianll2005l) . this reflection component introduces both 
line emission and a Compton hump to the spectrum. The emission 
line features are moderately strong (the Fe Ka line has equivalent 
width, EW ~ 140 eV for the upper spectrum in Figure[4]l and sig- 
nificantly broadened by the Doppler motions of the material (full- 
width at half-maximum, FWHM ~ 700 eV, determined by Gaus- 
sian fitting). In addition, Compton scattering of reflection line pho- 
tons in the outflow causes weak, red-skewed line wings to develop. 
This can be seen in the right-panels of Figure[5]which show details 
of the spectrum around the strong Fe Ka line. 

For moderate inclination angels (second row of Figure |5), 
the primary X-ray source is obscured by the high-ionization, fast 
outflow component. For these orientations, Compton scattering is 
the dominant opacity source along the line-of-sight although ab- 
sorption lines associated with high ionization state material (most 
clearly, Fe XXVI Lya) are also imprinted on the transmitted spec- 
trum. These features are generally rather narrow but blueshifted 
from their rest energies due to the net outflow along the line 
of sight. Emission features, typically slightly broader (FWHM 
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Figure 5. Spectra computed for three observer orientations (specified by /i ; 



9 where 8 is the inclination of the observer measured from the z-axis) 



for timestep 800 of the PK04 simulation. The left panels show the photon energy range 0.2 - 200 keV while the right panels show the 2-10 keV region 
(including the Fe K« line) in detail. In the left panels, the total spectrum (black) is subdivided into contributions from photons which have come directly from 
the primary X-ray source (red) and photons which arise from scattering/reprocessing in the outflow (blue). All spectra are normalized to the input primary 
X-ray spectrum (which was a pure power-law with photon index T = 2.1). The range of rest-energies for Fe, Ni and S K« transitions (6.4 - 6.97 keV, 7.5 - 
8.1 keV and 2.3 - 2.6 keV, respectively) are indicated in the upper right panel. 



~ 800 eV) and with somewhat more developed red wings than pre- 
dicted for polar inclinations, are also moderately strong for these 
orientations. 

At higher inclination, the absorbing column density is larger 
and the direct component of the spectrum becomes increasingly 
dominated by the scattered/reprocessed component (third panel of 
Figure[5}. This results in a spectrum with prominent emission lines 
(Fe Ka EW ~ 400 eV) and a stronger Compton hump. 

In most respects, the spectra we obtain from the two snapshots 
lead to qualitatively similar conclusions: in particular, the same 
characteristic spectra features (Compton hump, Fe Ka emission 
line, blended emission lines at soft energy, blue-shifted absorption 
lines) manifest and have similar dependence on the observer incli- 
nation. For fixed orientation, the emission features also have com- 
parable strength in the two snapshots, being just slightly stronger 
(typically by ~ 30 per cent in flux) in the second snapshot (re- 
call that the second snapshot has typically higher column densities 



which makes scattered/reprocessed radiation more dominant). In 
detail, however, there are some interesting quantitative difference 
between the absorption line properties of the two snapshots, mostly 
for intermediate observer orientations, which we will now discuss 
in greater detail. 



3.3.2 Complexity at intermediate orientations 

In both snapshots, the dependence of the spectrum on inclination 
angle is relatively modest when the orientation is either high or low. 
The sensitivity is strongest for intermediate orientations, mostly be- 
cause of the rapid change in the transmission properties of the wind 
for inclinations that pass through the complex outflowing structure. 
It is in this range of orientation angles that the most complex spec- 
tra occur but also where the least ambiguous outflow signatures 
(e.g. narrow, blueshifted absorption lines - see Paper I) can man- 
ifest. To map out the spectra in this sensitive regime, we repeated 
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the Monte Carlo simulations extracting spectra for a finely spaced 
grid of twenty inclination angles (with A9 = 1 deg) centred around 
the orientation used in the second panel of Figure [5] (48 deg < 9 < 
67 deg). A subset of the spectra obtained for this set of inclination 
angles in the timestep 800 snapshot is shown in Figure [6] divided 
into the component of radiation reaching the observer directly from 
the X-ray source and that formed by scattering/reprocessing in the 
disk wind structure. 

As expected, the direct component of the spectrum (upper 
panel in Figure[6]l is a stronger function of inclination than the scat- 
tered/reprocessed component. This can be readily understood since 
the scattered/reprocessed spectrum is formed by an integration over 
all positions in the flow and is therefore relatively insensitive to 
a small (~ 1 deg) change in inclination. In contrast, the attenua- 
tion of the direct component is determined only by the conditions 
in the narrow column of gas which obscures the small (assumed 
r e m = 6r 9 ) X-ray emission region. Since the properties of the 
line-of-sight column change rapidly with inclination, this compo- 
nent of the spectrum is strongly dependent on 8. Our calculations 
for the direct component of t he radiation field are in qualitative 
agreement with the results of ISchurch et alj d2009l) : they demon- 
strated that strongly orientation-dependent attenuation is expected 
in the PK04 simulation for 50 deg < 9 < 67 deg (see e.g. their fig- 
ure lfl But since their calculations effectively include only the di- 
rect component of the radiation field, we argue that their study over- 
estimates the role of absorption in shaping the X-ray spectrum: in 
our emergent spectra (composed of the the sum of the two compo- 
nents shown in Figure|6j, scattered/reprocessed light dominates for 
many inclinations leading to spectra with less dramatic absorption 
features and overall weaker dependence on the observer orienta- 
tion. 

The full spectra (direct plus scattered) around the Fe Ka re- 
gion for intermediate inclination angles are shown for both of our 
snapshots in Figure [7] This illustrates the most striking difference 
between the spectra obtained from the two snapshots, namely the 
behaviour of the blueshifted absorption lines. Owing to the fairly 
narrow range of polar angles at which a fast wind is present in the 
model, it is only for a modest fraction of lines-of-sight that clean, 
sharp absorption lines appear in our computed spectra for either 
snapshot. In the timestep 800 snapshot, an Fe Kq absorption line 
is present for a range of inclinations A9 ~ 12 deg. For the later 
snapshot, the range is even smaller, A9 ~ 5 deg. In our previous 
work (Papers I and II), narrow absorption lines also only appeared 
for only a minority of inclinations although they were generally 
more common than found here. This is likely a consequence of the 
simplicity of the velocity law adopted in the parametrized models 
(Papers I and II) - by adopting a smooth outflow velocity at all 
points, those models increase the fraction of lines-of-sight which 
pass through moderately opaque material with a significant com- 
ponent of velocity directed towards the observer. In contrast, the 
model considered here has a much more complex velocity field 
leading to a rarer occurrence of clean outflow signatures. 

Most importantly, the blueshifted absorption features are dra- 
matically different between the two snapshots. In the timestep 
800 snapshot, the Fe Ka absorption line not only manifests for a 

1 Note that perfect agreement between our calculation of the transmitted 
spectrum and those of Schur ch et al 1 120091) is not expected since our cal- 
culation includes the effects of ionization by scattered light and we have 
adopted a different spectral shape for the source of ionizing radiation. In 
general, we obtain a higher ionization state of the absorbing gas for inter- 
mediate inclination angles than they found. 



wider range of inclinations but it is generally significantly sharper, 
deeper and less blueshifted compared to the same inclination in 
the timestep 955 snapshot. The maximum Fe Ka blueshift for the 
first snapshot corresponds to only ~ 0.015c while it is as large as 
~ 0.06c for the later snapshot. The largest Fe Ka absorption EW is 
similar in both snapshots ~ 70 eV, although this occurs for different 
^-values (0 ~ 66 deg for timestep 800 and 9 ~ 60 deg for timestep 
955). For fixed inclination angle, the absorption EW changes quite 
significantly (a factor of two or more) between the two snapshots 
for most orientations. Thus the model predicts that the properties 
of blueshifted absorption features should not only be a strong func- 
tion of observer orientation but will also be significantly variable 
on timescales comparable to the time difference between our snap- 
shots (At ~ 5 years). We note that the calculated continuum level 
is also different between the two snapshots, most obviously around 
9 < 65 deg. This is attributable to the larger column densities for 
these lines-of-sight in the timestep 955 snapshot (see Figure[2j. 

In contrast, the model predicts that emission features (in par- 
ticular, Fe Ka) should be present in the spectrum for all observer 
orientations and that their character will be less dramatically time- 
variable (except for 9 > 75 deg, the Fe Ka emission flux typically 
changes by no more than ~ 30 per cent between the two snapshots). 
The relative insensitivity of the emission line flux arises from the 
fact that they are formed over an extended region in the flow and 
are thus less affected by details of the structure along the observer's 
line of sight. As noted above, it is at the same inclination angles for 
which absorption line features form that the emission lines are most 
intrinsically broad (FWHM ~ 0.8 keV) and also where they de- 
velop the most noticeably red-skewed wings via Compton scatter- 
ing in the flow (see Figure[8j. We note that the red line wings found 
in the current simulations are somewhat less well-developed than 
in the simplified models considered in Papers I and II; this is ex- 
pected since the fast wind component, in which Doppler shifts can 
most effectively give rise to the red-skewed wing, occupies only a 
relatively small region of the model. 



4 SUMMARY, DISCUSSION AND CONCLUSIONS 

We have applied our Monte Carlo radiative transfer code (see 
Papers I and II) to compute X-ray spectra for snapshots from 
a radiation-hydrodynamics simulation of a line-driven AGN disk 
wind (PK04). In most important respects, we find that the results 
of these simulations support conclusions drawn from our previ- 
ous studies of parametrized outflow models. In particular, (1) we 
find that a disk wind should imprint a wide range of spectroscopic 
features in the X-ray band and (2) the role of the wind in reflect- 
ing/scattering/reprocessing radiation in these simulations is at least 
as important as the part it plays in imprinting absorption signatures. 

Since we have considered only two snapshots from one model, 
one should not expect that the simulated spectra will quantitatively 
match observations. Rather the synthetic spectral features should 
be interpreted as broadly indicative of those which a disk wind 
may cause. This should guide future studies in which the hydro- 
dynamical modelling will be improved and the system parameters 
explored. It should also be borne in mind that the PK04 simulation 
was primarily an attempt to build insight to the problem of line- 
driven AGN winds and is not expected to capture all the detailed 
physics. In particular, since the model is axisymmetric, it will not 
correctly describe any small scale clumping in the flow. 

In the following sub-sections we first discuss the relevance of 
our results to the interpretation of X-ray spectra (Section f4. It . We 
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Figure 6. Direct (upper) and scattered/reprocessed (lower) components of the spectrum computed for 48 deg < 9 < 67 deg for timestep 800 of the PK04 
simulation (each inclination is indicated by a different colour; see colour bar in figure). The complete model spectrum (see Figure|7] left) is composed of the 
sum of these two components. Note the different ordinate scale used in the two panels. All spectra are normalized to the input primary X-ray spectrum (which 
was a pure power-law with photon index T = 2.1). 



then comment on the comparison of our spectra to those obtained 
from the parametrized models used in Papers I and II (Section[4j2} 
before discussing the implications of our work for the modelling of 
line-driven winds (Section [4.3b and, finally, highlighting some of 
the outstanding questions in Section l4~4l 

4.1 Implications for the interpretation of X-ray spectra 

4.1.1 Absorption lines 

Our first application of a Monte C arlo code to the modelling of the 
X-ray spectra of AGN Jslmll2005h was motivated by the identifica- 
tion of narrow, bl ueshifted absorption lines in the Fe K band (e.g. 
iPounds et al]|2003l) . Such features are the least ambiguous signa- 
ture of a fast, highly-ionized outflow. This study confirms that such 
features can arise in a plausible disk wind geometry. Our calcu- 
lations suggest that narrow, blueshifted Ka absorption lines will 
manifest in systems with line-driven disk winds for a limited range 
of orientations (very roughly Af? ~ 2 — 10 deg, centred around 
moderate inclination angles of 9 ~ 57 deg). For an isotropic dis- 
tribution of observed orientations, this would suggest that only a 
modest fraction of sources (somewhere around ~ 3 — 15 per cent) 



should show any such absorption lines. The fraction in the observed 
sample of X-ray bright AGN would be expected to be somewhat 
higher since very high inclination systems are likely excluded - 
but it should still be a minority. Thus the absence of this "smoking 
gun" signature of outflow is not evidence against the presence of a 
fast wind and the calibration factor required to convert the fraction 
of objects with observed narrow line features to the fraction with 
powerful disk winds may be large. 

Moreover, comparing the spectra obtained from the two 
timesteps we considered demonstrates the absorption line proper- 
ties should be significantly variable on timescales comparable to 
At ~ 5 years. Thus, even if the flow is persistent, re-observations 
of particular objects are not expected to yield very similar absorp- 
tion line properties over timescales of years: for a fixed observer 
orientation, the line blueshifts, widths and equivalent-widths can 
all change significantly. This is consistent with observations of e.g. 
PG121 1+143, an object where blueshifted Fe Ka absorption has 
been claimed but in which the absorption properties have varied 
between multipl e observations made over the course of the last 
decad e (see e.g. lPounds et al.ll2003l.r2005l ; |Pounds & Reevesll2007l 
2009). Our current radiative transfer simulations do not constrain 
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Figure 7. Spectra in the Fe Kq region for 48 deg < 6 < 67 deg (each inclination is indicated by a different colour) for the two snapshots (timestep 800 and 
955; left and right respectively), which are separated by At ~ 5 years. All spectra are normalized to the input primary X-ray spectrum (which was a pure 
power-law with photon index T = 2.1. The small-scale fluctuations are due to Monte Carlo noise in the simulations. 
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Figure 8. The spectrum of escaping Fe Kq photons computed for an ob- 
server inclination around 55 degrees for the timestep 800 snapshot from 
the PK04 simulation. The solid black line shows the complete spectrum 
of escaping Monte Carlo packets whose last physical interaction, excluding 
Compton scattering, was with an Fe Ka transition. The coloured lines show 
the subdivision of this profile into contributions from packets which under- 
went particular numbers of Compton scattering events after their interaction 
with an Fe Kq transition. The spectra are normalized to the peak of the line 
and have been slightly smoothed to suppress Monte Carlo noise. 



the level of absorption line variability expected on much shorter 
timescales (e.g. on the scale of typical X-ray exposures of several 
kiloseconds). Detailed consideration of shorter timescales will re- 
quire a more extensive, tightly-spaced sequence of snapshots and 



investigation of the role played by both small-scale structure in the 
flow (which will not be well-represented in the axisymmetric simu- 
lation considered here) and intrinsic, short time-scale variability of 
the primary X-ray source. We would like to note, however, that (i) 
variations in the wind structure itself occur on a very wide range of 
timescales, at least as short as ~ 10 5 s - thus we expect that some 
degree of variability will also be present on much shorter timescales 
than the time difference between the two snapshots used here; and 
(ii) short time-scale variability (timescales < 10 5 s) is much more 
likely to manifest in the absorption lines than any emission fea- 
tures since the light-crossing time for the simulation (~ 0.01 yr 
= 3 x 10 5 s) sets a lower limit for the timescale on which the 
emission from the outflow should show significant variability. 



4.1.2 Emission features 

Whilst narrow absorption features appear for only a limited range 
of observer orientations, other spectroscopic signatures associated 
with emission/reflection by the wind are predicted to appear for all 
inclination angles: 

(i) A Compton hump is present for all observer inclinations. 
This generally causes the spectrum to smoothly rise across the en- 
ergy interval ~ 10 — 25 keV. At harder energies it either flattens 
or, for high inclination angles, turns down. In the current simula- 
tion, the amplitude of the rise is quite modest, typically several tens 
per cent (measured relative to a power-law extrapolated from lower 
energies). 

(ii) For all observer orientations, the 2-10 keV band shows a 
fairly broad Fe Ka emission line (and corresponding, but weaker, 
Ni and S features). This feature is shaped by fluorescence, recombi- 
nation and scattering around the surface of the dense "failed wind" 
region and in the fast outflow. It can be fairly strong (emission EWs 
of at least ~ 150 eV for most orientations) and can develop a mod- 
est red-skewed wing due to Compton scattering of line photons in 
the outflow. Thus, for a realistic disk wind configuration, a signif- 
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icant component of strong, broad Fe Ka emission with a complex 
line profile can be easily produced. 

(iii) Below ~ 1 keV an additional emission component forms 
from multiply blended line and continuum features. These are dom- 
inated by reflection/reprocessing by the abundant elements Fe, O 
and C and a modest component of bremsstrahlung emission from 
the wind. Most of the emission features are strongly blended al- 
though there is always a fairly prominent O emission feature at 
0.7 keV and significant Fe emission at 0.8 keV. We note that, be- 
cause it is often a complex blend, the exact strength and shape of 
this low energy emission will be sensitive to changes in the relative 
element abundances (which were fixed to their solar ratios here). 

Qualitatively, each of the three classes of feature mentioned 
above is known from observations of AGN (hard band excesses, Fe 
Kg e mission, excess soft band emission; see e.g. lTurner & Milled 
2009). Given that there are many simplifying assumptions in the 
PK04 simulation and that we have made no attempt to tune the 
simulation parameters, we find it encouraging that such observable 
signatures are immediately predicted. 

Although the EWs of our Ka emission lines are comparable 
to those measured for the strongest narrow Ka line components 
in AGN, our em ission lines are too broad to account for these: 
IShu et al.1 J20ld) typically find Fe Ka line cores that have FWHM 
< 10, 000 km s" 1 while our lines have FWHM ~ 30, 000 km s -1 . 
The observed narrow line cores are therefore more likely to origi- 
nate in structures further out than the main line-forming regions in 
our simulations. Our synthetic emission line profiles can be more 
realistically connected with moderatel y broadened line c omponents 
such as reported in MCG-5-23-15 bv lBraito et aIU2007h . For such 
a case, the agreement in FWHM of the profile is much closer al- 
though the emission EWs obtained from the current model are too 
large by a factor of a few. Although interesting, the amplitude of 
the rise in flux between ~ 10 and 25 keV is small compared to 
that required to ex plain Suzaku obser vations of hard excesses in 
AGN (e.g. PDS 45 6lReeves et~all 2009 ; NCG 405 1 iTerashima etal] 
120091 : 1H04 19-577 iTurner et all 2009). To explain such strong ex- 
cesses requires absorption by high column-density optically-thick 
clumps along the line of sight. To fully explore this possibility in 
disk wind models will require full three-dimensional simulations to 
be undertaken in the future (see Section l4!4l . 

We argue that the qualitative properties of the spectra obtained 
from the current model (with no tuning/exploration of parameter 
space) encourages further study of disk wind models. We are cer- 
tainly forced to the strong conclusion that a line-driven disk wind 
should imprint several X-ray spectroscopic features in addition to 
"smoking gun" blueshifted absorption line features. Thus, if the fast 
outflows required to explain blueshifted absorption lines are associ- 
ated with a disk wind, it seems inevitable that the outflow will also 
have an important part in shaping many other spectral features, even 
when it does not directly obscure the primary X-ray source region. 

4.2 Comparison to parametrized models 

All the classes of spectral feature described above appeared in our 
calculations with simply-parametrized outflow models (see figs. 
3 and 4 of Paper II). Qualitatively, the features have comparable 
strengths and fairly similar shapes to those obtained with Model B 
in Paper II. The narrow absorption line features also have approxi- 
mately the same maximum blueshifts. The most apparent quantita- 
tive difference is the relative paucity of narrow absorption features 
in the spectra computed here. Also, despite having a wind mass-loss 



rate (Al w ) which is a significantly smaller fraction of the Eddington 
accretion rate, the model presented here predicts emission features 
which are comparably strong to those in the parametrized models 
used in Paper II (M w ~ 0.06A/Edd here but was > O.BAfEdd in the 
models adopted in Paper II). Thus the effective emission measure of 
the hydrodynamical wind must be higher than that of the equivalent 
parametrized model with the same mass-loss rate. In addition, the 
red-skewed emission line wings are typically less well-developed 
than in the parametrized models. 

These differences can all be attributed to the greater complex- 
ity of the hydrodynamical model: in contrast to the parametrized 
models, it is not a steady-state flow but has density and veloc- 
ity variations that are not simple functions of azimuth and radius. 
Moreover, neither the "failed wind" region nor the dense mate- 
rial around the disk at large radii are captured in the parametrized 
model. Both these regions contribute to the Fe Ka line emission 
(see Figure|4j i.e. these regions increase the effective emission mea- 
sure of the wind). 

We conclude, therefore, that the simply-parameterized mod- 
els provide a good means by which to study the influence of a disk 
wind on the spectrum. In particular, such models provide a reason- 
able description of how a successfully launched outflowing compo- 
nent imprints blueshifted lines and associated broad emission fea- 
tures from highly ionized material. There is, however, an important 
and obvious caveat with parametric modelling, namely that it is 
only useful to the degree that it reflects the physical properties of 
the wind. In our case, the parametric model is a poor representation 
of the "failed wind" component, a component which does have an 
affect on the overall spectrum. One must always be cautious about 
quantitative results derived from a parametric wind description. 

Taken together, our studies (Papers I, II and this work) suggest 
that models in which a smoothly outflowing component is domi- 
nant will predict relatively common blueshifted absorption features 
(and more strongly-developed red-skewed line wings) compared to 
a flow with complex topology. In particular, narrow blueshifted Ka 
absorption is most likely to be clearly seen in flows with moderate 
covering fractions along lines-of-sight that have modest continuum 
optical depths (see discussion in Paper I). Relative to emission fea- 
tures, these absorption line features are likely to be rarer in sce- 
narios in which other flow components (e.g. the "failed wind" in 
the simulation presented here) have significant covering fractions. 
Thus, there are good prospects that an unbiased census of X-ray ab- 
sorption line features could, in the future, place constraints on the 
typical geometry required for wind models. 

4.3 Implications for the theory of AGN winds 

In this study, we have adopted the existing radiation- 
hydrodynamics simulation of PK04 and performed detailed 
radiative transfer simulations as a post-processing step, as- 
suming that the underlying density/velocity distributions are 
appropriate. Our detailed X-ray radiative transfer calculations, 
however, point to a potentially important shortcoming of the 
radiation-hydrodynamics simulation. 

In the radiation-hydrodynamics calculation it was assumed 
that X-ray radiative transfer in the wind was a purely absorptive 
processes. This leads to the presence of an extended wedge around 
the base of the wind which is only weakly ionized since it is 
shielded from direct irradiation by a centrally located source. It is 
in this low ionization material that the radiative line force due to the 
uv disk photons is most effective in accelerating the flow. However, 
we have explicitly shown that scattered/reprocessed radiation has a 
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pivotal role in determining not only the emergent X-ray spectra but 
also the ionization/thermal conditions in a disk wind. In particular, 
with our multi-dimensional simulations we find that much of the 
low density region behind the shielding gas is still rather signifi- 
cantly ionized thanks to scattered light (see Figure^. 

Additional ionization can reduce the line driving, pote ntially 
leadin g to a weaker wind (or even quenching the wind; see IProgal 
2005). We speculate that the likely consequence of accounting for 
ionization by scattered X-rays would be that a larger fraction of 
the material would be in the "failed wind" component whilst less 
material will be successfully accelerated to escape velocity. To in- 
vestigate this quantitatively will require further developments of 
our radiation-hydrodynamics code and additional simulations. We 
note, however, that for sources accreting close to Mmi, the forma- 
tion of the "failed wind" component seems unavoidable and our 
results suggest that even this alone is an important structure for the 
formation of the X-ray spectrum. Thus it seems unlikely that struc- 
tures formed due to the pressure exerted by the uv radiation field of 
the disk do not play any part in formation of the X-ray spectrum of 
AGN accreting close to the Eddington limit. 



4.4 Further work 

Important future improvements to the modelling should include in- 
vestigating how the outflowing gas might continue to affect the 
spectrum at larger radii and lifting the assumption of axi-symmetry. 
A fully three-dimensional flow could be more complex and would 
likely include more density inhomogene ities (see e.g. right pan- 
els of fig. 2 in iKurosawa & Progal l2009). This may influence the 
observables. In particular, optically thick cl umps might have a 
role in explaining hard-band excesses (see e.g. lReeves et al.ll2009l ; 
iTurner et alj|2009h . Since both the wind structure and primary X- 
ray source luminosity are likely to be time-variable, the radiation 
transport should ultimately be make time-dependent to allow for 
the calculation of time-lags between spectral features forming in 
different regions and departures from thermal/ionization equilib- 
rium. 

Here, our study has been limited to line-driven winds 
(as simulated by PK04) but other processes may also con- 
tribute to the launching of AGN outflows: e.g. magnetically 
driven or thermal winds from either an irr adiated disk (e.g. 
Luketic et af] l2O10h or an AGN torus (e.g. iDorodnitsvn et alj 



2008al lbl: IDorodnitsvn & Kallmanl |2009h . Thus, in the future, 
we will extend our studies to consider alternative wind sce- 
narios and investigate whether there are prospects to discrimi- 
nate between possible launching mechanisms based on spectro- 
scopic observations. In particular, it will be important to com- 
pute spectra based on a magnetically driven disk wind model. 
It has been shown that magnetic fields could drive a disk 
wind over a relatively wide range of radii and could explain 
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